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Abstract — In cellular networks, admission control and beam- 
forming optimization are intertwined problems. While beam- 
forming optimization aims at satisfying users' quality-of-service 
(QoS) requirements or improving the QoS levels, admission 
control looks at how a subset of users should be selected so that 
the beamforming optimization problem can yield a reasonable 
solution in terms of the QoS levels provided. However, in order 
to simplify the design, the two problems are usually seen as 
separate problems. This paper considers joint admission control 
and beamforming (JACoB) under a coordinated multicell MISO 
downlink scenario. We formulate JACoB as a user number 
maximization problem, where selected users are guaranteed to 
receive the QoS levels they requested. The formulated problem is 
combinatorial and hard, and we derive a convex approximation 
to the problem. A merit of our convex approximation formulation 
is that it can be easily decomposed for per-base-station decen- 
tralized optimization, namely, via block coordinate decent. The 
efficacy of the proposed decentralized method is demonstrated 
by simulation results. 

Index Terms — admission control, distributed optimization, 
downlink beamforming 

I. Introduction 

Coordinated beamforming (CoBF) |[T1. |[2l is a recently 
studied technique to mitigate intercell interference (ICI) in 
the downlink of multicell cooperative systems. In CoBF, the 
neighboring BSs share the same frequency band and employ 
beamforming for data transmission. The transmit beamformers 
at different BSs are coordinately designed according to the 
channel conditions and certain design formulation, e.g., maxi- 
mum system throughput, minimum transmit power, to name a 
few. Compared to fully multicell cooperative techniques such 
as network MIMO HI, CoBF has an advantage that the BS 
cooperation overheads are not as significant, and yet appealing 
performance may be achieved. 

Meanwhile, admission control also plays an important role 
in cellular systems. As cellular systems are usually congested, 
with lots of users awaiting service, it is necessary for the BSs 
to decide which user is served or not. Admission control refers 
to methods of selecting users. While admission control and 
beamforming are commonly seen as two separate problems, 
they are fundamentally dependent on each other. Recent work 
has demonstrated that by considering admission control and 
beamforming jointly, promising system performance can be 
achieved O, H. However, joint admission control and beam- 
forming (JACoB) is a challenging problem. It has been shown 
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that JACoB is NP-hard even under a single-cell scenario 0. 
Hence, as a compromise, one may consider approximation 
approaches. 

This paper describes a JACoB approach to CoBF in multi- 
cell MISO downlink. A distinguishing part of the present work 
is that our proposed JACoB formulation can be easily decom- 
posed for decentralized implementation, and the decentralized 
process is considered even more straightforward than those in 
CoBF (without admission control), e.g. J2], J5j- A key idea of 
our approach is to use the now popularized l\ approximation 
method. We very recently note that in a concurrent work H, 
the authors have studied t\ approximation for joint admission 
control and power control (i.e., no beamforming). Neverthe- 
less, the work |4| does not investigate the multicell CoBF 
scenario and, more importantly, decentralized optimization 
considered here. It is also worthwhile to mention Q, which 
considers JACoB under a single cell scenario. The formulation 
used there is based on a mixed-integer program formulation, 
and is processed by semidefinite relaxation (SDR). While it 
is not difficult to see that the idea of (3) can be extended 
to the multicell scenario, one needs to assume centralized 
optimization and presently there is no reported work on how 
the method in [3| can be decentralized. In our simulations, we 
will show that our decentralized method yields a performance 
quite on a par with the centralized method in (3). 

II. System model 

Consider a cellular system with M coordinating BSs. Each 
BS is equipped with N transmit antennas. In each cell, there 
are K single-antenna user terminals; thus the total number 
of users in the system is KM. The set of users associated 
with the zth cell, or the ith BS, is denoted by tCi C K. = 
{1,2, ...,KM}. We assume that K = Kill- •■U/C M and /Cjfl 
ICj = whenever i ^ j, i.e., each user is served only by one 
BS. The scenario of interest is downlink, with an emphasis on 
CoBF. Assuming that the BS-to-user channels are frequency- 
flat and slow, and that the linear unicast transmit beamforming 
scheme is employed, we can characterize the CoBF system 
performance by the received signal-to-interference-and-noise 
ratios (SINRs) (for more complete system model descriptions, 
see the literature, such as J2)): 



SINR„=- 



melC i{q) \{q} 



H 
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(1) 



where q £ {1,2,..., KM} is the user index, i(q) £ 
{l,...,M} denotes the BS with which the qth user is asso- 
ciated (i.e., i(q) is such that q £ /Q( g )), hj, q £ C N is the 
channel response from the jth BS to the qth user, a q is the 
noise variance, and w q £ C w is the beamforming vector for 
the qth user. 

The difference between the SINR model in (HJ and a single- 
cell-based SINR model is that the former explicitly models the 
ICI, which is given by the third term in the denominator of 
((T). On the contrary, in the single-cell case, the ICI is usually 
treated as a constant, and is absorbed by a q . The idea of CoBF 
is therefore to ask the BSs to coordinately design {w q }^(, 
so that ICI may be jointly mitigated. 

To motivate the study of JACoB, let us quickly review a 
CoBF design problem, specifically, the design proposed in 0. 
In that design, the BSs jointly design the beamforming vectors 
{ w q]^=l suc h that the SINR of each user is no less than a 
user-requested threshold j q . Moreover, the design objective 
aims at minimizing the total transmit power. This amounts to 
the following design optimization problem: 

q =i Klla 

s.t. SINR, >7 9 , q = 1,2,..., KM, (Z) 

where the last constraints in (O are per-BS power budget 
constraints, with P m ax.i specifying the maximum transmit 
power of the ith BSLil and || • ||2 is the £2 norm. At first 
glance problem © seems to be nonconvex, since the SINR 
constraints are nonconvex in {w q } q J( . Actually, problem (f2]i 
can be solved in a convex and tractable fashion, using ei- 
ther the second-order cone programming formulation or the 
semidefinite relaxation (SDR) formulation; see J2) and (6) for 
more detail. 

Cellular systems are usually congested, with lots of users 
awaiting service. A subsequent issue relevant to the CoBF 
problem (0 is that we may be unable to find a beamforming 
solution {w q } q Jf that satisfies all the users' SINR requests. 
In other words, problem (f2]i may be infeasible. To illustrate this 
issue, we simulated the feasibility rate of problem (f2]i against 
the total number of users KM. The simulation result is plotted 
in Figure Q] The feasibility rate was evaluated by counting the 
number of instances for which (f2]i is feasible, under randomly 
generated channels. We observe that problem (f2]i has a low 
feasibility rate when the number of users is large. 

III. Joint admission control and beamforming 

This work considers joint admission control and beamform- 
ing (JACoB). The problem is stated as follows: 

Joint admission control and beamforming (JACoB): 

Select a maximum number of users, such that there exists 
a beamforming solution {w q } that satisfies all the selected 
users' SINR requests. 



'As a minor point to note, the previous study [2] does not incorporate the 
per-BS power budget constraints. 
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Fig. 1: The feasibility rate of the CoBF problem (0. M = 3, 

N = 8, 7, = 6 dB, and P max ,i = 46 dBm. 

In the following, we will first provide an optimization 
formulation for JACoB, and derive a convex approximation to 
the formulated problem. Then, a decentralized method based 
on the convex approximation will be developed. 

A. Centralized method for JACoB 

Our endeavor starts with formulating JACoB in a mathe- 
matically convenient form. First of all, let 

W q =w q w q 1 , q=l,...,KM, 

and observe the following equivalence (cf. (Q])): 

SINR, > lq <=► > mWm}^), (3) 

where we define 

/,(w™) = i+ E H H ^( £ Wm )) 

+Tr(iJ % ), 9 ( E W m -±W q )), 

me/C iW \{g} lq 

and Hj lQ = hj q h^ q / 'a q . We claim that JACoB can be 
formulated as the following £0 minimization problem: 

min ||t||o + e£fi/Tr(Wg (4a) 

s -t- E geKs TrW) < P max ,i, i - 1, M, (4b) 

W q h 0, V q, (4c) 

t q = max{0, /,({W m }™)}, V q, (4d) 

rank(W g ) = 1, V q, (4e) 

where < e < 1/ YhU Pmax, 

i is a penalty parameter, W q y 
means that W q is positive semidefinite, and ||t||o is the £0 
norm, which counts the number of nonzero elements in t. 

Let us describe why problem (|4]i delivers the above defined 
JACoB goal. Firstly, constraints ( f4cb and ( |4el ) are equivalent 
to W q = w q w^ . Secondly, by substituting ( |4d| i into the first 
term of (l4al i. i.e., ||t||o, and observing (J3J, we can see that 
||t||o is counting the number of unserved or unadmitted users. 
Hence, if we ignore the second term of (f4ab . then problem (01 
minimizes the number of unadmitted users. Thirdly, the second 
term of (|4ab , i.e, e ^2^ =1 Tr(W q ), is a penalty term. It is used 
to encourage more power-efficient beamforming solutions. It 
can be shown that problem (|4]i with e = (i.e., direct 



unadmitted user minimization) achieves the same number of 
unadmitted users as problem (0 with e < l/5^ i=1 Pmax.i- 

Problem (0 is difficult to solve. As a remedy, we adopt a 
convex approximation approach. Our approximation involves 
two steps. First, we replace the hard £q norm function by the l\ 
norm, which is now a popularized trick in compressive sens- 
ing. Second, we remove the rank-one constraints d4el ). which 
is well known as SDR J6). The above two approximations lead 
us to the following l\ approximate JACoB problem: 



min ||*||i+e£f-f TrW) 

s - t - Y,q&Ki Tr ( W q) ^ Pmax,i, * 
t, > 0, t, > /,({W m }^l), 
^ 9 b 0, V q, 



1,...,M, 



(5a) 

(5b) 
(5c) 
(5d) 



where || ■ ||i is the t\ norm. Note that in (0, we replace d4dl > 
by d5cT ). which can be easily verified to be equivalent. 

The l\ approximate JACoB problem <[5j is convex. In fact, 
problem (0 can be written as an SDR Hence, for centralized 
implementation, we can solve problem (0 by using a readily 
available SDP solver. Moreover, we show that the second 
approximation, i.e., SDR, is a tight relaxation: 

Proposition 1 For e > 0, any optimal solution {W*}^ 1 / of 
problem (0 must satisfy rank(W*) < 1 for all q. 

The proof of Proposition 1 is skipped here owing to the limit 
of space. The idea behind the proof is to examine the KKT 
conditions and exploit the rank-one structure of Hi( q v q . We 
should also note that Proposition 1 is different from the SDR 
tightness results in (71, which may be seemingly similar at 
first look. Simply speaking, Q studies more general Hj q 
(which may take any rank), but may not solve the problem in 
Proposition 1. Proposition 1 means that solving problem © 
automatically leads to a set of beamforming solutions (recall 
W q = WqW^ for rank-one positive semidefinite W g ), and 
there is no loss in applying SDR. 

B. Decentralized method for JACoB 

A significant advantage of the £% approximate JACoB 
formulation in (0 is that it can be easily decomposed for 
decentralized optimization. To see this, let 

m = {W m } meKt , i = l,-,M. 

Notice that Wi corresponds to the beamforming vectors 
controlled by the ith BS. Now, by substituting (Bel into (Bat , 
we can reformulate (0 as 



KM 

min > 

tw \ KM 

iyV «'9=l q=l 



max{0, f q (Wi, Wm)} + e.Tr{Wq] 



s-t- MWq) < P„ 



i = l,...,M, 



(6) 



WqhO, q= 1,...,KM. 

A unique feature with problem (0 is that the constraints are 
per-BS decoupled (note that this is not the case with the 
CoBF problem in (0). As a result, we can directly apply 



per-BS alternating optimization. To be specific, we employ 
block coordinate descent (BCD). In BCD, we update only one 
beamforming block Wi, while holding the other blocks fixed. 
This BCD update is done cyclically with respect to the BSs, 
until some stopping rule is satisfied. 

A curious question is whether the above-described BCD 
method would converge to the optimum of problem (0. Unfor- 
tunately, this may not be guaranteed — BCD may not converge 
to the optimum for problems whose objective functions are not 
continuously differentiable, even if the problem is convex J8). 
The function max{0, x} seen in problem (0 exactly falls into 
this case. To remedy this, we apply a smooth approximation 
to (0 using the one-sided Huber function 

{0 , if x < 0, 

0.5a; 2 , if < x < 1, 
x - 0.5 , if x > 1. 

The Huber function h(x) is continuously differentiable in x. 
Applying the approximation max{0, x} w h(x), we obtain the 
Huber approximate JACoB problem: 

KM 

mi $ M E (M/ 9 (V^,...,WA/))+eTr(Wg) 



9=1 



S-t. X/ijgACi Tr(Wq) < Pmax,ii 

W q t 0, q= I,..., KM. 



i = 1,...,M, 



(7) 



From this point on, we will concentrate on the BCD of 
problem (0. 

Let us consider the BCD update of problem (0 with respect 
to the ith block Wi, holding the other blocks {Wjjj^i = 
{^mjmttj fixed. The respective problem is 

KM 

, w m r E M/«(w\, {Wjh^i)) + E ^r(w m ) 

1 q—1 mGfCi 
S -t- Y^meKi Tr (^m) < Prnax,i, W m h 0, TO £ £j, 

which can be expressed as a convex problem (see, e.g., J9)): 

KM 

mm > -ul + v„ ) + V eTr(W m ) (8a) 



q—1 m£/Ci 

s.t. u q + v q > /,(W„ {W,}^), V q, (8b) 

Uq, Vq>0,\/ q, (8C) 

E ^ P max ,i, W m h 0, m e Ki. (8d) 

Constraints d8bl seems to indicate that full knowledge of 
{VVj jj-ij is required, in order to solve the BCD update (0. 
Actually, this may be not necessary. Notice that for q G JCi, 
constraint d8bl can be expressed as: 



Ug + Vq 



>i+E^+ Tr (^( E w - 



j¥* meK,\{q} 

and for q ^ JCi, 

Uq + Vq > 1 + E + T^-^i,? ( E W " 



7g 



where Ctj, q are scalar constants defined as: 



Algorithm 1 Block coordinate descent method for 



Hence, if the ith BS knows i) the matrices {Hi iq } q Jf , i.e., the 
channel response from the ith BS to the users in the system; 
and ii) the scalar constants {&j,q}jj&i,q£ic, then problem (0 
can be solved independently at the ith BS. In fact, the first 
premise can be satisfied automatically as it can be assumed 
that each BS knows the channel response from itself to the 
users in the system 10. To satisfy the second premise, we 
can utilize the backhaul link between the BSs. Specifically, 
the scalars {&i !q } q£ )c can be computed and broadcast to the 
other BSs after the ith BCD update is solved. There are KM 
real numbers to be broadcast at each iteration. This justifies 
our claim that the ith BS can solve (0 alone. The BCD method 
for is summarized in Algorithm [TJ 

Recall that the reason for employing the Huber function in 
is to provide a smooth approximation to JACoB, avoiding 
the original nondifferentiable objective function which may 
result in BCD non-convergence problems. But can the smooth 
approximation guarantee convergence to the optimum? By 
invoking an available BCD convergence analysis result IflOl , 
we have the following claim: 

Fact 1 ifTOl The sequence {{W g }k generated by 

Algorithm [7] has limit points and every limit point of the 
sequence {{W q }k is an optimal solution to 0. 

Readers are referred to Proposition 6 in iTTOl for more detail. 
We should note that the important premises for us to use this 
available result are that the objective function of is convex, 
continuously differentiable in {Wi}^i and the constraint 
set for each Wj is convex and compact. Furthermore, by 
extending Proposition 1, we can prove that the intermediate 
solutions {W q }q£tCi in Algorithm 1 at the ith BCD update 
are always of rank-one: 

Proposition 2 For e > and for each i, any optimal solution 
{Wq}q^c t of the ith BCD update in Algorithm 1 must 
satisfy rank(Wg ) < 1 for all q in /Q. 

IV. Deflation heuristic 

Both l\ and Huber approximate JACoB problems (cf. prob- 
lems and (01, respectively) can be seen as some kind of 
"soft decision" formulations for handling admission control. 
In order to select more users for service, we can apply a hard 
decision using the deflation heuristic. Similar to J3), |4), the 
heuristic is initialized by considering all users in the system, 
then the users are dropped one-by-one. At first the BSs solve 
(0 or either centrally or using the BCD method. Our user 
dropping rule is based on the value of t* = max{0, f q (-)} 



Require: initialization - {W q } q eK- 
l: k = 1; 

2: For each i = 1,2,..., M, the ith BS computes {Cli tq } q& ic 

and broadcasts them to the other BSs; 
3: repeat 

4: for i = 1 to M do 

5: The ith BS solves ® given {Qj,q}jjH,qEic to obtain 
{Wq }q£Ki- The scalars {£li jq } qe ic are computed 
and broadcast to the other BSs; 

6: end for 

7: k = k + 1; 

8: until convergence. 

9: return an optimal solution to - {W q k ^} qe /c- 



which relates directly to the satisfiability of the SINR threshold 
for user q. The user with the largest t* will be dropped. 

When the number of users in the system is too large, we 
may encounter cases where the optimal solution to or is 
trivial, i.e., W q = for all q. Here we state an easy-to-check 
condition for identifying such cases. 
Fact 2 (Prescreening condition^ If 

* 9 ({%), m W))^o,v?elC ! (9) 

then solving or gives a trivial solution, i.e., W q = 
for all q in JC, where 

*g({-Hi(g),m}me/c)) — d + E ^i(.9),m ^*(«),«- 

The proof is omitted due to space limitation. Inspired by 
fact 2, we now adopt a prescreening procedure where we 
drop the users gradually until condition (0 gets violated. 
Specifically, at each time, we remove user q with & q (-) that 
gives the largest minimum eigenvalue As condition (0 can 
be checked in closed-form, the prescreening procedure can be 
run at a low complexity. The deflation heuristic, together with 
prescreening, are summarized as followfll: 



2 A similar condition has been discovered recently in [4:| for the joint power 
and admission control problem. Our results applies to the case with CoBF. 

3 Note that both the deflation heuristic and prescreening procedure can be 
operated in a decentralized manner. 



Deflation heuristic: 

• Initialize: a set of users requesting service - K, = 
{1,2,..., KM}. 

1) (Prescreening) Check condition (0. If it holds, 
then remove user m from /C according to m = 
argmaXq A TO j„(*q(-)) and repeat 1). Otherwise, go 
to 2). 

2) (Deflation) Solve (0 or for {i* W*} q£K . If t* = 
for all q <E JC, terminate. Otherwise remove user m 
from JC according to m = argmax g t*, repeat 2). 

• Return: a set of selected users JC and rank-one ma- 
trices {W q }q£ic that decomposes into beamforming 
vectors satisfying the SINR requirements. 
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Fig. 2: Total no. of users admitted, M = 3, N = 8, A' = 15. 

V. Numerical results 

This section presents numerical results for the proposed 
JACoB methods. The simulation environment is similar to 0. 
We focus on a multicell scenario with 3 coordinating BSs 
where each BS is separated from the others by 2.8 km. For 
each simulation trial, the users' positions and their respective 
channels are randomly generated. The users are separated from 
their respective BS by at least 0.7 km and are assigned to the 
nearest BS. The channel is assumed to experience both small- 
scale and large-scale fading. The noise variance a 2 q is -92 dBm 
and the receive antenna gain is 5 dBi. There are 8 transmit 
antennas and 15 users are assigned to each BS (i.e., 45 users 
awaiting service in total). 

Here, the centralized method and the decentralized method 
refer to the deflation heuristic which uses a centralized solver 
for (0 and the decentralized BCD method for (0, respectively. 
The scalar e is chosen as 10~ 5 and the BCD method terminates 
when the relative change in objective value is less than 10~ 2 . 

Figure|2]shows the performance of different JACoB methods 
in terms of the total number of users selected for service. Here, 
two benchmarking methods are compared. The "SDR deflation 
method" is adopted from [3 | and modified to operate in the 
multicell scenario. Note that it is a centralized method which 
may not be decomposed straightforwardly. The "single cell 
deflation method" is a modified algorithm from [11], which is 
originally proposed as a suboptimal alternative to the CoBF 
problem (UJ. The main feature of ifTTI is that the ICI levels are 
always constrained below a fixed threshold, and therefore the 
beamforming design can be done independently at each BS. 
The subsequent single-cell beamforming design is processed 
by our JACoB method. Note that the resultant method can be 
implemented in a per-BS decentralized manner without any 
BS coordination. 

Turning back to Figure [2] we compare the performance 
of with the proposed centralized method. We observe 
that there are some performance gains with the proposed 
method in the high SINR regime. Furthermore, the proposed 
decentralized method achieves a performance on a par with 



TABLE I: Total no. of iterations, M = 3, N = 8, A = 15. 



Threshold (7) 


With prescreening 


Without prescreening 


12 dB 


53.080 


61.790 


20 dB 


31.030 


59.290 



its centralized counterpart. The decentralized method should 
also be compared to the single cell deflation method, where 
the numerical results has clearly demonstrated the benefits of 
allowing BSs coordination. 

Table Q] demonstrates the efficacies of the proposed decen- 
tralized BCD method with prescreening procedure in terms of 
the total number of iterations. The total number of iterations 
is defined as the total number of BCD iterations consumed 
throughout the deflation heuristic, where multiple instances 
of (01 are solved. Note that the load on the backhaul link is 
directly proportional to the iteration count. The iteration counts 
reported in Table U confirms that significant reduction in the 
number of iterations can be achieved with the prescreening 
procedure. 

VI. Conclusion 

The contributions of this paper are twofold. First, we have 
developed a formulation of joint admission control and beam- 
forming (JACoB) for coordinated multicell downlink, wherein 
an efficient convex approach is proposed. Second, we have 
built a decentralized JACoB method via a simple BCD pro- 
cedure. Simulation results have shown that the decentralized 
method achieves a performance on a par with the centralized 
method with fast convergence. 
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